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In this paper we study kinetically rough surfaces which display anomalous scaling in their local 
properties such as roughness, or height-height correlation function. By studying the power spectrum 
of the surface and its relation to the height-height correlation, we distinguish two independent causes 
for anomalous scaling. One is super-roughening (global roughness exponent larger than or equal to 
one) , even if the spectrum behaves non anomalously. Another cause is what we term an intrinsically 
anomalous spectrum, in whose scaling an independent exponent exists, which induces different 
scaling properties for small and large length scales (that is, the surface is not self-affine). In this 
case, the surface does not need to be super-rough in order to display anomalous scaling. In both 
cases, we show how to extract the independent exponents and scaling relations from the correlation 
functions, and we illustrate our analysis with two exactly solvable examples. One is the simplest 
linear equation for molecular beam epitaxy , well known to display anomalous scaling due to super- 
roughening. The second example is a random diffusion equation, which features anomalous scaling 
independent of the value of the global roughness exponent below or above one. 

PACS numbers: 5.40.+j, 05.70.Ln, 68.35.Fx 
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I. INTRODUCTION 

The problem of the dynamics of growing surfaces has 
attracted very much attention in recent years. A large 
variety of discrete models and continuous equations have 
been investigated numerically and analytically in order to 
explain the diversity of rough surfaces observed in non- 
equilibrium natural processes In many cases, it is 
observed that the surface morphology displays time and 
space fluctuations which have scaling properties akin to 
those of equilibrium critical phenomena. Most of the 
work is thus devoted to identifying the different univer- 
sality classes to which the models studied belong, which 
are basically determined by the exponents characterizing 
the fluctuations of the interface. 

We will be interested in this work in the dynamics of a 
one dimensional surface, which will be described by the 
stochastic variable h(x, t) giving the height of the surface 
at time t above substrate position x. The physical prop- 
erties of the system are thus described by the moments 
of h(x,t). For instance, the global width 



where the scaling function f(u) behaves as 



W(L 1 t) = (((h(x,t)-h(t)) 2 ) 1 J 2 ), 



(1) 



is a measure of the interface fluctuations. In ([!]), (...} x 
denotes spatial average over the whole system of lateral 
size L, h(t) = (h(x, t)) x is the average value of the height 
at time t, and (...) is an average over realizations of the 
noise. Starting from a initially flat surface, the global 
width is observed to satisfy in many cases the dynamic 
scaling Ansatz of Family- Vicsek S, 



/(«) 



const, if u ^> 1 
u x if u C 1 



(3) 



The roughness exponent x characterizes the surface mor- 
phology in the stationary regime, in which the correla- 
tion length £(t) ~ t 1 ^ (z is the so called dynamic expo- 
nent) has reached a value larger than the system size L. 
This happens for times larger than the saturation time, 
t ^> t s (L), which scales with L as t s (L) ~ L z . The ratio 
j3 = x/ z is called growth exponent and characterizes the 
short time behavior of the surface. 

The existence of dynamic scaling as in (0) , plus the as- 
sumption of the self-affine character of the interface — in 
the sense that there is no characteristic length scale in the 
surface besides the system size, and therefore all scales 
obey the same physics — allows one to obtain the expo- 
nents from the local width w(l,t), which measures the 
surface height fluctuations over a window of size I <C L. 
At times larger than t s (l) ~ l z the local width is thus 
expected to saturate and 



(4) 



w(L,t) = t^ z f(L/t 1 / z ), 



(2) 



where the local roughness exponent xioc should equal the 
exponent obtained in (||J3|), xioc = X- F° r short times, 
t <C t s (l), the local width scales with time as w(l, t) ~ t 13 , 
which is the short scale analog of the scaling behavior for 
the global width given in (||). It is an interesting fact to 
stress here that the equality xioc = X m ((§) is n °t guaran- 
teed in general when Family- Vicsek scaling holds for the 



1 



global width (§), since the self-affinity of the interface is 
an additional independent condition. 

Although the above scaling picture has been success- 
fully applied to a great variety of models, which includes 
for instance the importat universality classes associated 
to the Edwards- Wilkinson (EW) Q and Kardar-Parisi- 
Zhang (KPZ) equations, it is not always valid. In 
very recent studies of discrete models and continuum 
equations, mainly related with the growth of surfaces by 
molecular beam epitaxy (MBE), it has been found that 
in the case of the so-called super-rough surfaces, i.e. for 
surfaces with a global roughness exponent \ > 1, the 
usual assumption of the equivalence between the global 
and local descriptions of the surface is not valid. In these 
systems the local width (or equivalently the height- 
height correlation function, see below) does not saturate 
as expected in the standard Family- Vicsek scaling, but 
crosses over to a new behavior in the intermediate time 
regime l z C ( « L z , characterized by a different growth 
exponent /3*, w(l,t ^> l z ) ~ I t 13 * , where /3* = (3 — 1/z. 
It is easy to see that now the local width saturates only 
when the global width does (this is, t s (l) = L z for all I) 
and w(l,t ^> L z ) ~ IL X_1 in saturation, which, accord- 
ing to (Q), yields a local exponent xioc = 1 for super- 
roughening. This scaling behavior which is not encom- 
passed by the the Family- Vicsek Ansatz has been termed 
anomalous scaling in the literature [p|JTo|] and later differ- 
ent models have been studied in which xioc and /3* take 
values different from 1 and (3 — 1/z, respectively jLl[ . 

While it is currently clearly understood that super- 
roughening always leads to anomalous scaling, in a recent 
paper jl2| two of us have demonstrated that growth mod- 
els in which % < 1 may also exhibit an unconventional 
scaling behavior with similar scaling relations between 
exponents. In that reference, for an analytically solv- 
able growth model with tunable values of x, we identified 
anomalous dynamic scaling with the lack of self-affinity. 
Also a class of discrete MBE growth models studied nu- 
merically in Ref. [jllj behaves anomalously, although the 
interfaces generated were not super-rough at all. That is, 
surfaces can display anomalous scaling no matter what 
their value of the global roughness exponent is, either 
X > 1 or x < 1- bi the presence of anomalous scaling, 
not all length scales are equivalent in the system; the 
scaling behavior is different for short (local) and large 
length scales, hence xioc and x differ. 

Our aim in the present paper is to complete the pic- 
ture presented in Ref. JlJ] , showing that the mechanisms 
which lead to anomalous scaling behavior can be sepa- 
rated into two classes, according to the behavior of the 
structure factor or power spectrum of the surface, S(k, £), 
to be defined below. One of the mechanisms is super- 
roughening, even if S(k,t) behaves as in the Family- 
Vicsek case. The other mechanism corresponds to what 
we term intrinsic anomalous scaling of the structure fac- 
tor, where an independent exponent appears which mea- 
sures the difference between the short and large length 
scale power laws, namely, the difference between the local 



and global roughness exponents xioc and x- We will show 
how to identify the corresponding anomalous scaling by 
extracting the independent critical exponents from the 
correlation functions. As stated above, for simplicity we 
restric our further analysis to growth models in one di- 
mension, but the conclusions are expected to be valid in 
arbitrary dimensions. We note that, although the scaling 
relations to appear below for the super-rough and intrin- 
sic anomalous cases exist already in the literature [jnjU > 
seemingly the independence between the occurrence of 
anomalous scaling and the value (larger or smaller than 
one) of the global roughness exponent x has not been 
adequately understood so far. Since a clear picture of 
the different scaling behaviors which can appear is cru- 
cial to correctly identify them in experiments and models, 
we believe a somewhat detailed presentation is in order. 
It is attempted in Section II. The rest of the paper is 
organized as follows. In Section III we illustrate the con- 
clusions of Sec. II by an analysis of two growth models 
formulated in terms of stochastic equations. Even though 
the two models are exactly solvable, we choose to study 
them numerically in order to illustrate how our scaling 
analysis can be applied to more general cases. The first 
model examined is the linear MBE equation, known to 
display anomalous scaling due to super-roughening. The 
second model we consider in Sec. Ill is the stochastic dif- 
fusion equation with random diffusion coefficient studied 
in fl3| , [l2|| , which exemplifies a case of intrinsic anom- 
malous scaling. Finally, we present some conclusions in 
Sec. IV, and discuss consequences for scaling analysis of 
rough surfaces in experiments and numerical simulations. 
Appendix A contains details on the calculations discussed 
in Sec. II. 



II. ANOMALOUS DYNAMIC SCALING 

In order to study a problem of kinetic roughening the 
height-height correlation function 

G(l,t) = (((h(l + x,t)-h(x,t)) 2 ) x ), (5) 

is frequently used. This correlation function scales in 
the same way as the square of the local width, G(l,t) ~ 
w 2 {l, t), and provides an alternative method to determine 
the critical exponents. 

However, as we will see in the following, a more trans- 
parent understanding of the complete dynamic scaling 
can be obtained by studying the Fourier transform of the 
interface height in a system of linear size L (see e. g. Jj[] 
and references therein), 

h(k,t) = L- 1/2 J2[h{x 1 t)-h{t)}exp{ikx), (6) 

X 

where the spatial average of the height has been sub- 
tracted. In this representation, the properties of the sur- 
face can be investigated by calculating the structure fac- 
tor or power spectrum 
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S{k,t) = (h(k,t)h(-k,t)), (7) 

which contains the same information on the system as 
the height-height correlation function G(l,t) defined in 
(||), both of them being related by Jl4j 

G ( l ^) = J [1 - cos(fci)]S(M) (8) 

27T / L<k<TT I 'a 
1""/°- dh 

oc / — [1 -cos(kl)]S(k,t). 

J2it/L 27r 



L — > 00 having been taken. Analogously one can arrive 
at the local scaling given by 

G(!,t) ~ { % /Z :| Jvl J J } = l 2x gFvil/t 1 ^), (13) 

where the scaling function gFv{u) is constant for u <C 1 
and gFv(u) ~ w _2x for u»l. This is the usual situa- 
tion in which ^ < 1 and the local and global roughness 
exponents are equal (x = Xioc and /?* =0). 

-Z. Super-roughening 



A. Family- Vicsek scaling 

Family- Vicsek scaling reads, when expressed in terms 
of the structure factor, 



S(k,t) =k-^ +1 h FV {kt l ' z ), 
with sfv the following scaling function 



SFv{u) 



const. 



u 



2X+1 



if u > 1 
if u < 1 



(9) 



(10) 



Indeed, Eqs. (^,|l^) can easily be seen to be equivalent 
to Eqs.(||, ||), by noting that the global width is nothing 
but the integral of S(k, t), i.e. 

W 2 (L, t) = ~ £ 5 ( fc < *) = / |^ fc < *)' C 11 ) 

k 

where the momentum integral is limited to 2tt/L < k < 
n/a and represents a continuum aproximation to the sum 
over the discrete set of modes. In a discrete growth 
model, a is identified with the lattice spacing. Note as 
well that Eq. (^TJ) implies that, for kt 1 ^ 1, the spec- 
trum does not depend on time, and hence at saturation 
(t 1 /* > L), S(k,t) is a pure power law independent of 
system size. 

Going back to real space, and having assumed a 
Family- Vicsek scaling for S(k,t), Eq. (|l0|), one obtains, 
using Eq. (|), 

G(l,t) = £[l-cos(W)] SFV J k x t! Z) dk 

= l 2 * [; [1 - cos (u)]-^0±du, (12) 



where we have introduced the ratio p = t x / z /I. The lim- 
iting behaviors of (|T^) are very simple to study when 
X < 1; the details can be found in the Appendix. For 
instance, at saturation t x l z L, the scaling function in 
the integrand of Eq. ( |l2] ) can be substituted for its (con- 
stant) limiting behavior at large arguments (see ([To|)), 
and one is left with G(l,t) ~ l 2x , the limits a — > 0, and 



In the case of growth models generating super-rough 
surfaces (x > 1), but with a structure factor fulfilling 
Family- Vicsek scaling, the integrals in ( [l2| ) are divergent 
in the limit I <C t x / z (p 1) for L — » 00, given the 
strong singularity at the origin of integration j^] . Taking 
the limit I <C t x l z first for fixed a, L (see Appendix), one 
obtains 



G(l,t) 



I 2 ^(x- 1 )/ 2 if I < t 1 ^ < i 

Z 2 L 2(x-l) if / « i < tl/* 



(14) 



so that x; oc = 1 and f3* = — 1/z. In the early time 
regime t 1/z < / G(Z,i) — t 2x/z . As has been re- 

marked in Refs. 1 15 I3] , the fact that xioc cannot exceed 
1 for super-rough surfaces (x > 1) is a purely geometric 
property which follows from definition (^|) . 



B. Intrinsic anomalous roughening 

Apart from the super-roughening case, there is another 
important situation which leads to anomalous scaling. In 
some growth models — an example of which is discussed 
in Sec. Ill — , the structure factor may present an uncon- 
ventional scaling not described by Family- Vicsek (|lC|). 
Let us consider a dynamic scaling form for S(k,t) as in 
Eq. (^J) but with the scaling function s fv (u) replaced by 



s A (u) 



,20 



if u > 1 

i 2 * +1 if U < 1 



(15) 



where the label A denotes intrinsic anomalous spectrum. 
Here 9 is a new exponent which "measures" the anomaly 
in the spectrum. In a system of size L, Eqs. (||), ( |l5| ) 
hold only up to the saturation time t s (L) = L z , after 
which the system size L replaces the correlation length 
t l l z in all expressions. Thus in particular, at saturation 
the structure factor depends on the size of the system 
as S(k, t) ~ L B k e ~( 2x+1 \ As a consequence, the station- 
ary spectrum shows severe finite size effects, to the extent 
that it is not defined in the thermodinamic limit L — > 00. 

A scaling behavior such as Eqs. @, ((TJ) for the struc- 
ture factor does not affect the behavior of the global 
width, which preserves its Family- Vicsek form, W(L,t) ~ 
for t < L z and W(L, t > L z ) ~ L x . On the contrary, 
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the local properties of the surface change dramatically if 
S(k,t) scales as in (|lj]). That can be seen by computing 
the height-height correlation function from Eq. (g), 



G(l,t) = l 2x 



, si S A {pu) 

cos [u)\ — , , du 



(16) 



which as before (and see Appendix) gives G(l,t) ~ t 2/3 
for times t <C l z ■ However, for intermediate times 
I 2 « ( < L z (/>> 1) the integral ( p6| ) picks up a non- 
trivial contribution from the behavior of sa(u) at large 
arguments, so that 



G(Z,t) 



j2 X 28 = px-W-flS/z 



(17) 



Thus, the complete scaling of the height difference corre- 
lation function (or, equivalently, the square of the local 
width) can be written as 



G(l,t) 




if t 1 '* « I « L 
if I < t x l z < L 
if !<£<tV« 



ptgAQ/t 1 '"), 



(18) 



where the local roughness exponent is Xioc = X — 0, 
[3* = 9/z = j3 — Xloc/ z j an d the scaling function <?a(w) is 
noi constant anymore for small arguments, but behaves 
as 



9a(u) 



u -2( X -Xloc) if U< 1 

2 * if u»l 



a 



(19) 



The fact that 9 ^ in Eq. (y_5[) yields a local roughness 
exponent xioc ^ X an d an anomalous growth exponent 
/3* 7^ 0. Therefore, there are now three independent ex- 
ponents describing the scaling properties of the surface, 
whereas for Family- Vicsek scaling (even in the presence 
of super-roughening) there are only two. Which triplet of 
independent exponents one considers depends mainly on 
the kind of scaling properties one is measuring. In exper- 
iments normally it is local properties of the surface (local 
width, etc.) that are measured, and one might naturally 
consider as independent exponents e.g. Xloc, z, and /?». 
On the other hand, the focus of numerical simulations 
are usually finite size effects, thus probing global prop- 
erties of the surface such as the global width or power 
spectrum. In the latter cases a common choice is e.g. 
X, 0, and (3. We should stress that, although the scal- 
ing relations for exponents are exactly the same as for 
the super-roughening case discussed above, the present 
type of anomalous scaling is due to 9 ^ and may and 
does occur in growth models with % < 1. The scaling 
behavior (|l^) of the local width is formally equivalent to 
that of super-roughening (for which case xioc = 1 and 
(3„ = [3 — l/z), a fact which has produced some confu- 
sion in previous works where both anomalies have been 
identified in some way. 

In view of the local scaling, Eq. (jl^), the structure 
factor can be rewritten conveniently as 



, tPx+i)/* if fctV*«i 

MMJ~<i k -(2 Xloc + l) t 2(x-X^)/z if fctV*»l ' ^> 



where we can observe two interesting facts that charac- 
terize what we call an intrinsic anomalous scaling. First, 
S(k,t 1 ^ z 3> fc -1 ) decays as k~^ 2xioa+1 \ and not follow- 
ing the fc - ( 2x+1 ) law characteristic of Family- Vicsek scal- 
ing (super-roughening case included). Second, there is 
an unconventional dependence of S{k,t 1 / Z ^> k^ 1 ) on 
time which leads to a non stationary structure factor for 
t x l z k ^> 1. The combination of these two facts allows 
to distinguish between anomalous scaling due to super- 
roughening (x > 1) an( i intrisic anomalous scaling. 



III. MODELS AND NUMERICAL RESULTS 



In this section we present numerical simulations of 
two different growth equations in order to illustrate the 
theory we have just discussed in the previous section. 
First, we discuss the linear MBE equation as a paradig- 
matic example of Family- Vicsek behavior of the struc- 
ture factor, in the presence of super-roughening, lead- 
ing to anomalous scaling of the height-height correlation 
function. Second, we study a random diffusion model in 
which one can tune the values of the global exponents 
X, /?• We compare the exact solution of the model with 
numerical integrations of the equation of motion. By an- 
alyzing the structure factor, we will see that the second 
model constitutes an excellent example of what we have 
termed intrinsic anomalous scaling in the previous sec- 
tion, independent of the value of x below or above one. 



A. Linear MBE equation: super-roughening 



Let us start studying the simplest equation relevant for 
MBE growth Jl8|,pp[|. For a one dimensional interface it 
reads 



dh _ d 4 h 



rj{x,t), 



(21) 



where 77 is a white noise with zero mean, and the constant 
K > 0. Exponents can be easily found by dimensional 
analysis, = 3/8, z = 4, and x = 3/2, thus Eq. @ 
generates super-rough surfaces. This equation has been 
widely studied by many authors; however for the sake 
of clarity we wish to sumarize here briefly some numeri- 
cal results for it that illustrate our scaling analysis in an 
example of super-roughening. 

We have solved numerically Eq. ( |2l| ) by means of 
a finite differences scheme. To study the structure 
factor, we have calculated S(k,t) in systems of L = 
16,32,64,128,256. 
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FIG. 1. The inset shows the structure factor for the linear 
MBE equation, Eq.(21), calculated in a system of total size 
L = 128 for times t = 10 2 ,3 x 10 2 ,5 x 10 2 ,7 x 10 2 ,9 x 10 2 . 
The solid straight line has slope —4 and is ploted to guide the 
eye. In the main panel the data are collapsed using exponents 
X = 3/2, z = 4. The solid line has slope 4, showing that the 
scaling function of the structure factor is consistent with Eq. 
(10). 



An inspection of the inset of Fig. 2 shows that the local 
width on length scales I <C L does not saturate at times 
t 3> l z but continues to grow in time as w(l, t ^> l z ) ~ t@* 
with /3* ~ 0.13, in agreement with the theoretical predic- 
tion @ for Eq. @ /3„ = p - 1/z = 0.125. 

As discussed above, in the presence of anomalous scal- 
ing the height-height correlation function (or the local 
width) satisfies a dynamic scaling form, although with 
a scaling function which differs from the Family- Vicsek 
one, see Eqs. (HI), ©, and ©. As stressed in 0] 
(for the height-height correlation function), the plot of 
w{l,t)/l x vs. ljt x l z allows one to determine the form of 
the scaling function since 

w{l,t) 



g{i/t 1,z ). 



(22) 



The data in the inset of Fig. 2 have been collapsed for 
1 = 5,..., 150, and plotted in Fig. 2. The scaling func- 
tion g(u) obtained is not constant for m«1, but scales 
with its argument as u~( x ~ Xloc \ with xioc. = 1 for this 
case. This is consistent with the scaling behaviour Jl^), 
or equivalently with (|18|) for xioc = 1- 



In the inset of Fig. 1 numerical results are shown 
for L = 128, and 200 realizations of the noise. For 
k^ 1 3> t x l z , the structure factor decays as k~ 4 (i. e. 
the global roughness exponent \ = 3/2) and does not 
display any anomaly. In Fig. 1 we plot a collapse of the 
data for the exponents z = 4 and x — 3/2. As we see, 
the scaling function of the structure factor indeed has the 
Family- Vicsek shape (|H)]) |l^] . Thus there is no anomaly 
in the behavior of S(k, t) for Eq. (pl|). However, there is 
anomalous behavior if we look at the local properties of 
the surface in real space. We have calculated the local 
width, w(l,t) for different window sizes I in a system of 
total size L = 1000, and averaged over 15 realizations of 
the noise. 




0.0 ,. 

log 10 (l/t 1 ") 

FIG. 2. Numerical results for the local width of the linear 
MBE equation in a system of size L = 1000. In the inset we 
plot w(l,t) vs. time for two different window sizes I = 5 (o) 
and I = 150 (□). The straight lines have slopes 0.13 (solid) 
and 0.37 (dashed), corresponding to (5* and respectively. 
Data are collapsed in the main panel for I = 5, up to / = 150 
using exponents % — 3/2, z = 4. The straight lines correspond 
to u~ x (dashed) and u ~^ x ~ Xloc ^ (solid) and have slopes —3/2 
and —1/2 respectively. 



B. Random diffusion model: Intrinsic anomalous 
scaling 

Next let us consider the more complex case of intrinsic 
anomalous scaling. We focus our attention on a growth 
model with a random diffusion coefficient governed by 
the equation 

^=JK Dw ^ t) ) +,M) ' (23) 

The diffusion coefficient D{x) > is distributed accord- 
ing to the probability density P(D) = N^D^^ f c (D/ D c ), 
where the parameter cj> characterizes the intensity of the 
disorder. is merely a normalization constant and the 
cutoff function is f c (y) = 1 for y < 1 and f c (y) = 
for y > 1. If (j> < 0, disorder in the diffusion coeffi- 
cient does not play any role and the Edwards- Wilkinson 
H exponents x = Xioc = 1/2, /3 = 1/4 are recovered. 
Thus the disorder is termed weak. On the contrary, for 
strong disorder, < <f> < 1, the critical exponents XiP 
are disorder-dependent (through the value of <p) . 

There are two major features that make this model in- 
teresting. On one hand, it is linear and can be solved 
exactly |l3|,[l2| . As a consequence, exponents are depen- 
dent on the parameter in a known way. In particu- 
lar, the global roughness exponent x takes a continuous 
range of values from x = 1/2 to x = oo when var ying 
the intensity of disorder from = to 1. In Ref. |12j | 
it has been found that the random diffusion model al- 
ways exhibits anomalous scaling for < <\> < 1. Hence 
anomalous scaling in this model occurs for surfaces with 
roughness exponent x > 1 well as for those with x < 1. 
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Therefore, this model constitutes an excellent example of 
anomalous scaling not due to super-roughening. 

As stated above, the analytical study of (|3|) allows to 
determine exactly the critical exponents (see jL2|). For 
the region of interest here, < < 1, one has 







1 



1 



2(2 



2(1-0)' 



Xlc 



(24) 



It is a curious property of model (|23J) that xioc remains 
constant when changing the intensity of the disorder, in 
contrast to the wide range of variation existing in \- This 
is a feature also found for many models of rough epitaxial 
growth pT[ |. In Sec. II we have demonstrated that in the 
case of anomalous kinetic roughening the time exponent 
/?» is always given by the scaling relation /?» = f3 — Xioc/ z - 
Thus, from (p4|), /3» depends on as /3* = 0/(4 — 2<f>) for 
a given disorder parameter 0. 

We have performed simulations of the discretized ver- 
sion of Eq. ( p3| ) in a system of lateral size L = 1000 for 
different degrees of disorder in the strong disorder phase, 
< < 1. In Fig. 3 we plot w(l,t) for two values of 
0. For = 2/3, the exact global roughness exponent is 
larger than one, \ = 3/2, and one observes in the upper 
panel of Fig. 3 that the scaling behaviour of the local 
width is anomalous with /3* = 0.27 ± 0.03 in agreement 
with the theoretical value /?* = 1/4. 




2.0 3.0 4.0 
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FIG. 3. Local width vs. time for the random diffusion 
model, Eq. (23), calculated over windows of size I = 5 (o) 
and I = 150 (□) for two different degrees of disorder. In the 
upper panel (0 = 2/3), the straight lines have slopes 0.27 
(solid) and 0.38 (dashed), to be compared with the exact val- 
ues /?» = 1/4 and p = 3/8, respectively. In the lower panel 
(0 — 1/3), the lines have slopes 0.12 (solid) and 0.31 (dashed), 
to be compared with /3* = 1/10 and /3 = 3/10, respectively. 
Thus, the scaling of model (23) is anomalous (i.e. 0* 7^ 0) 
independently of the value of the global roughenss exponent. 

In the case of a disorder parameter 0=1/3 (see lower 
panel of Fig. 3), for which % — 3/4 < 1, the scaling of 
w(l,t) is also anomalous with an anomalous time expo- 
nent /3* = 0.12 ± 0.04 to be compared to /3* = 1/10 pre- 
dicted by the exact analysis. Thus, in equation ( p3| ) the 
existence of anomalous scaling is completely independent 
of the surface being super-rough or not. 

As we did above for the linear MBE equation, we have 
determined the scaling function g(u) of the local width 



by the data collapse shown in Fig. 4, in which the ex- 
ponents obtained from Fig. 3 have been used. For both 
degrees of disorder = 1/3 and 2/3, the scaling function 
behaves as expected for anomalous scaling, and has an 
analogous shape to that found for the local width in the 
super-rough case. 
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FIG. 4. Data collapse of the results displayed in Fig. 3. In 
the upper panel (0 = 2/3), the exponents \ = 3/2, z = 4.08 
have been used. The straight lines are plotted as a guide 
to the eye having slopes —3/2 (dashed) and —1 (solid). For 
0=1/3 (lower panel), the data show a good collapse for 
X = 3/4, z — 2.42. The straight lines have slopes —0.25 
(solid) and —0.75 (dashed). In both panels the scaling agrees 
with formula (25). 



The complete scaling function g(u) defined in 
thus given by 

-(x-xioe) if u < 1 



9(u) 



if u > 1 



(25) 



where x is the disorder-dependent global roughness ex- 
ponent given in ( p4|) and xioc = 1/2. 

The intrisic anomalous character of the scaling in this 
model can be diagnosed through the use of the structure 
factor. We have calculated S(k,t) in systems of sizes 
L = 16, . . . , 512. Figure 5 shows our results for L = 128, 
and 200 realizations of the disorder. 
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FIG. 5. Structure factor of the random diffusion model 
Eq. (23) for two different degrees of the disorder and times 
t = 10 2 ,3 x 10 2 ,5 x 10 2 ,7 x 10 2 ,9 x 10 2 . Upper panel, re- 
sults for = 2/3 (x = 3/2). Lower panel, data for 0=1/3 
(x = 3/4). The shift in time reflects the intrinsic anomalous 
character of the scaling. 

For both values of the disorder parameter = 1/3, 2/3 
the spectrum decays as k~( 2xioc+1 ' ) (not as fc~( 2x+1 )) and 



G 



is clearly shifted for different times. This scaling be- 
haviour is the one in Eq. (^0|), which we associated in 
Sec. II. B with that of growth models having an intrisic 
anomaly. This can be better appreciated when collapsing 
the curves of Fig. 5 as shown in Fig. 6, which displays the 
S(k,t) data collapses for <fr = 1/3, 2/3, and yields a scal- 
ing function s(u) with a form consistent with Eq. (15), 
and not with the Family- Vicsek form, Eq. ([n|) jl9| . 




FIG. 6. Data collapse of the graphs in Fig. 5. The expo- 
nents used for the collapse in the upper panel (<j) — 2/3) are 
X = 3/2 and z = 4.08. The straight lines have slopes 2.2 
(solid) and 4 (dashed). In the lower panel (<j> = 1/3) the ex- 
ponents used for the collapse are \ — 3/4 and z = 2.42. The 
straight lines have slopes 0.6 (solid) and 2.5 (dashed). In both 
panels the scaling function is described by Eq. (15). 

As we see, this example clearly shows that intrinsic 
anomalous behaviour may occur in the case of a rough- 
ness exponent % = 3/2 > 1 (0 = 2/3) as well as for 
X = 3/4<1 (0-1/3). 



IV. CONCLUSIONS 

As we have seen in the previous sections, one of the 
main conclusions of the present work is that anomalous 
scaling of rough surfaces is a more general (and inde- 
pendent) phenomenon than those cases associated with 
super-roughening. In this sense several issues associated 
with anomalus scaling remain to be clarified. One of 
them is its relation and/or interplay with the phenom- 
ena of multifractality and intermitency, known to take 
place in the super-rough case |2^,^l| . Moreover, we have 
seen that the intrinsic anomalous surfaces display a novel 
type of spectrum thus far not considered for self-affine ge- 
ometries, which are usually defined as to have a Family- 
Vicsek spectrum p2| . An outstanding issue in this regard 
is to clarify the physical meaning of the exponent 8, as 
well as extending the currently available renormalization 
group techniques to be able to calculate the value of this 
exponent. On the phenomenological side, when perform- 
ing experiments or interpreting numerical simulations of 
discrete growth models, we believe the existence of in- 
trinsic anomalous scaling is a very relevant issue. On 
one hand, in the presence of intrinsic anomalous scaling 
the various correlation functions behave somewhat sim- 



ilarly to the Family- Vicsek case. However, indeed there 
are important differences. For instance, S(k,t) scales as 
/ c -(2xioc+i) ; an d not as k~( 2x+1 \ the graphs of w(l,t) and 
S(k,t) are shifted for different times, etc. There do exist 
data collapses, but again they are different from the usual 
Family- Vicsek type. For instance, if instead of plotting 
w(l,t)/l x vs. l/t 1 /* as recommended in jll] and done here 
in Figs. 2, 4, one chooses to represent w(l,t)/t^ vs. l/t 1 / 2 
as is frequently done, one gets, in the Family- Vicsek case, 



w(l,t) 



const, if l/t 1 / 2 < 1 
(ttW) X if l / tl/z » 1 



(26) 



whereas for intrinsic anomalous scaling one gets 



w(l,t) j const. if l/t 1 / 2 < 1 



(27) 



that is, in both cases there is data collapse with a very 
similar shape of the scaling function. However, the slopes 
of the corresponding scaling functions for large argu- 
ments arc different. In the Family- Vicsek case, the slope 
coincides with the value of the exponent assumed to 
achieve the collapse of the data. However, in the anoma- 
lous case it does not. Therefore, it is crucial to check 
whether the slope of the scaling function does or not coin- 
cide with the assumed exponents. We believe this has not 
always been done when analyzing data from experiments 
and/or numerical simulations, and may have added to 
certain confusion existing in the literature in the identi- 
fication of universality classes. Moreover, the intermedi- 
ate time regime existing in surfaces with anomalous scal- 
ing introduces difficulties in the evaluation of exponents 
through the common use of local measurements such as 
the local width, frequently employed in experiments and 
numerical simulations. This may lead to the assignment 
of erroneous effective values (see to the exponents 
of surfaces hypothesized to behave in the simple Family- 
Vicsek fashion, while a more accurate description of their 
true scaling properties may come through the use of the 
anomalous scaling Ansatz. 
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APPENDIX A: LOCAL SCALING PROPERTIES 
OF SURFACES 

In this Appendix we study the scaling properties of the 
height-height correlation function G(l,t) from the differ- 
ent behaviors of the structure factor S(k, t). In all cases, 
the starting point is the general formula 



G(l,t) 



[1 _ cos s { k f /Z J dl , 



fx 



[1 — cos (i 



s(pu) 



II 



2X+ 



j du, 



(Al) 



where s(-) is the corresponding scaling function of the 
structure factor under study, and p — t 1 '* jl. 



1. Saturation regime, i<Ij< t 1 ^ . 

In this regime, p 3> 1, a nd th e smallest possible value of 
the argument of s(pu) in ( |A1[ ) is also pnl/L ~ t 1 ^ /L 
1. Therefore we can replace the scaling function s(pu) 
by its behavior at large arguments, to get the following 
results: 

Family-Vicsek with x < 1 ; I n this case the scaling 
function s(pu) = Spv(pu) of Eq. (fol). We get 



G(Z,i) - P x 



1 



COST* 



,2x+l 



dw ~ l 2x , 



(A2) 



and 



where the limits L — > 00 and a — ► can be taken, 
the integral yields a finite constant. 

Family-Vicsek with x > 1: In this case the scaling func- 
tion is as above s{pu) = spvipu). However, for y > 1 
there is an ultraviolet divergence in the integral (A2). We 
can avoid the divergence by keeping finite the lower inte- 
gration limit (i.e. taking the limits a — ► 0, t 1 ^ /L — > 00 
for a fixed L). In this way we have 



G(l,t) 



I 2 * 
I 2 L 2 ^ 



1 



cos it 



du 



r 1 1 



2(1-X) 



(A3) 



where we have kept only the most singular term (as can 
be shown integrating by parts in the above formula). As 
stated in the text, this implies a local value for the rough- 
ness exponent xioc — 1- The case % = 1 is a marginal 
situation. Following the same steps as in Eq. (A3), one 
obtains 



G(l,t) ~ Z 2x logL, 



(A4) 



that is, xioc = X — 1 U P to logarithmic corrections. 

Intrinsic anomalous scaling: Now the scaling function 
we have to use in Eq. (Al) is s(pu) — sa(pu) ~ (pu) , 
see (£J|). Thus 



G(l,t) = l 2x p 



20 



1 



cosu 



u 



2 X +l-28 



du. 



(A5) 



The integrals are convergent as long as = x — Xioc > 
and xioc < 1 (the latter being always fulfilled as shown 
in 15 lq|). At saturation we thus have 



G(l,t)~l 2 ^L 2e = 1 2 *'°<L 



20 



(A6) 



hence the difference between local and global scaling be- 
havior (note, however, G(l ~ L,t) ~ L 2{ - e+ * l °^ = L 2x ). 



2. Early dynamics, t 1/z < I < L. 

In this case the m inim um value of the argument of the 
scaling function in (Al) is pirl/L ~ t x l z /L < 1, so that 
in principle the integral probes the whole scaling function 
s(pu), for both small and large values of its argument. In 
this subsection we first consider the case in which, fur- 
ther, we have p = t x l z jl <1. Let us perform a change 
of variables in expression ( |A1[ ) : 



G(l,t)=l 2 *p 2x 



TTt 1/z /a 



27Tt 1 /z/L L 



1 



y 

COS I - 



y 



s (y) , 



(A7) 



We want to study the asymptotics of the integral in (A7) 
in the limit p — > 0. Using standard techniques (see e. g. 
]l7t), and assuming lim a ^ s {y) < 00 j it is easy to show 
that, for p — > 0, 



Mldy + 0(p), (A8) 



pTTt 1 /* /a 




a: 


p CO 


1 — cos 










Jo 



where the integral appearing in the last equation is well 
behaved in its limits, and s(y) is any of the scaling func- 
tions we are considering, either spy(y) or s^(y). There- 
fore we get 

G(l, t) - l 2x p 2x + 0{p 2x+1 ) - t 2x ' z for t x ' z < /, 

(A9) 

that is, the scaling properties of the height-height corre- 
lation function in the early dynamics are the same irre- 
spective of the scaling behavior of the structure factor, 
whether Family-Vicsek or intrinsically anomalous. 



3. Intermediate time regime, I < t 1/2 < L. 

Finally, we study the case in which, analogously to the 
previous subsection, the minimum value for the argu- 
mentof the scaling function along the integration region 
in (Al) is small — so that the whole scaling function con- 
tributes in principle — but the parameter p = t L ' z jl is 
large. That is, we are in the intermediate time regime. 
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In this case, again we have to distinguish between FV 
and intrinsic anomalous scaling of S(k,t). 

Family-Vicsek sca ling with x < 1: When s(pu) = 
sfv(pu) i n Eq. flAl), the limits L — > oo, a — > can 
be taken after p 



(|Al|), the limits L 
— > oo, and 



G(f,t)~J^ / (1 



1 



I.2X+1 



- Z 2x , (A10) 



being the integral a finite constant. 

Family-Vicsek scaling with x > 1- Again for the case 
of super-roughening with FV s caling in S(k, t) an expan- 
sion such as that in Eq. QA10D leads to a divergent inte- 
gral. The behavior of G(l, t) is best analyzed in the form 
(A7). For x > 1j we can take the limits t 1 / 2 /a — > oo, 
f7*/L ^ to get 



G(i,i)~J a V*~ 



7/.VX+1 



/2 ^(x-l)/,_ (A11) 



In the marginal case x = 1, the integral in the above 
expression diverges when taking the mentioned limits 
t 1 / 2 1 a — > oo, t 1 / 2 /L — > 0, being the most divergent term 
logarithmic in time, that is: 

G(l,t) ~ / 2 logi. (A12) 

Comparing with Eq. this yields /3» = for x = lj 

meaning logarithmic corrections in time. 

Intrinsic an omalous scaling. When we use s(pu) = 
sa(pu) in (Al) and proceed as in (A1C), the behavior of 
the scaling function at large arguments now contributes 
with a leading term (in the limit p — > oo for fixed p l/a 
and p l/L) of the form 

r°° 1 
G{l,t) ~ Z 2 V 9 / (1 ~ cos M )^i- du ~ pix-V^/' 



hence the anomalous behavior of the height-height cor- 
relation function summarized in 
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